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Abstract 

We investigate the critical scaling behavior of finite systems in the canonical 
ensemble. The essential difference with the grand canonical ensemble, i.e., 
the constraint on the number of particles, is already known to lead to the 
Fisher renormalization phenomenon that modifies the thermal critical singu- 
larities. We show that, in observables that are not Fisher renormalized, it 
also leads to a finite-size effect governed by an exponent y\ that depends on 
the temperature exponent yt and the dimensionality d as y\ = — \2y t — d\. We 
verify this prediction by a Monte Carlo analysis of several two-dimensional 
lattice models in the percolation, the Ising and the 3-state Potts universality 
classes. 
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The central task of statistical physics is to calculate the partition function and the 
thermodynamic observables. Different ensembles may be employed to this purpose, which 
then naturally yield different forms of the thermodynamic functions. In the thermodynamic 
limit, the different approaches generally yield equivalent results for the relations between 
the thermodynamic variables. However, differences occur in finite systems. We illustrate 
this by means of a lattice gas, described by a finite size L, a temperature T and a third 
parameter governing the number of particles. In the canonical ensemble this third parameter 
is the particle density p. A calculation of an observable A in the canonical ensemble will 
thus, in principle, yield its expectation value (A) c as a function A^(T,p,L). In the grand 
ensemble, one employs the chemical potential p as the third parameter. One may thus 
calculate the expectation value (A) g = A^ 9 \T, p,L) as well as the grand canonical density 
(p) g = p^(T, p, L). In the thermodynamic limit one expects in general that 

A^(T,p,L) = A^(T,p^(T,p,L),L) (L -> oo), (I) 

for expectation values of the form (A) = A(T)P(T), where P(T) is the probability of 
state T in the pertinent ensemble. It does not apply to quantities obtained by differentiation 
of observables to T, such as the specific heat. This is evident when we define A' = dA/dT: 

A^(T,p,L) = A' {c \T,p^(T,p,L),L) + ^^. (2) 

The last term violates Eq. (I). It leads to the Fisher renormalization [1] effect that may 
affect even the exponents of leading critical singularities. 

In this paper, we point out that, in the canonical ensemble, new finite-size effects ap- 
pear also for quantities that do not involve differentiations to T or other temperature-like 
variables. They appear because substituting p^ in the right-hand side of Eq. (1) is not pre- 
cisely equivalent with taking the grand canonical expectation value. Namely, in the grand 
ensemble, p is still allowed to fluctuate, which is not the case in the canonical ensemble. 
We derive this new finite-size effect for a rf-dimensional lattice gas, described by a reduced 
Hamiltonian 7i with variables (Tj = (I) denoting the absence (presence) of a particle on 
lattice site i. The grand partition sum is 
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Z^(T,^L) = J2J2...J2eM-n], (3) 

(7l =0 <J2 =0 (Tjv =0 

where the sum is performed independently on all N = L d lattice gas variables. The particle 
density p follows from differentiation of Eq. (3): 

/m rx 1 dlnZfe) 

Other thermodynamic quantities A can be obtained similarly by differentiation to the conju- 
gate parameter. This leads to the following form for the expectation value of an observable 
in the grand ensemble 

AtofaT) = 5>(M)exp[-H]/Zfe). (5) 
M 

The sum on {a} is shorthand for the sums in Eq. (3). 
The canonical partition sum is 

z^( P ,T) = J2K^^M-n], (6) 

M 

where N p = Np is the particle number and the Kronecker S imposes the constraint. The 
canonical expectation value of A is 

A (c \p,T) = J2^ P ^ k A({a})eM-n]/Z^. (7) 

M 

Therefore, A^ and A^ are related as 

A y^ 1 )- z w Z ( g ) 

N 

= A^(N p /N,T)P(p,N p /N), (8) 

N p =0 

where the grand probability P(p,p) that a particle density p occurs is equal to the ratio 
Z^/Z^. In the thermodynamic limit, Eqs. (5) and (7) should be equivalent as long as p 
and p satisfy Eq. (4). 
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For simplicity, let A be a quantity whose leading singular finite-size-scaling term is a 
constant, such as a the Binder ratio [2] or other dimensionless finite-size amplitude ratios. 
Near a critical point p = p c , T = T c , its finite-size scaling behavior reads [1,3] 

A^(p, T c , L) = + J2a k (p- Pc ) k L^ + bV>\ (9) 

k 

where the last term allows for a finite size effect at T c . Eq. (8) thus becomes 

A®(iJL,T e ,L) =A^ + Y,ak((p-Pc) k )L^ + bLy\ (10) 

k 

where ((p — p c ) n ) = J2n p =o(p ~ Pc) n P{p-, N p /N) is the n-th moment of (p — p c ) in the grand 
ensemble. The second moment scales as 

((p-p c ) 2 )=r L- d + r 1 L^- 2d . (11) 

At criticality, the term with k — 1 in Eq. (10) is suppressed, and substitution of Eq. (11) 
yields 

A (g) (/i c , T c , L) = 4 c) + wL-\ 2yt - d \ + bL yi + ■■■ (12) 

where we have used y p = d — y t for 2y t — d > and y p = y t for 2y t — d < according to 
the Fisher renormalization [1,3]. The asymptotic value is A^ = A^f 1 for 2y t — d < and 
A^f = A^f 1 + a 2 ri for 2y t — d> 0. Eq. (12) should match the finite-size scaling formula of A 
in the grand canonical ensemble. Here we are on the firm ground of the scaling properties 
of the free energy, in terms of scaling fields associated with the intensive thermodynamic 
parameters. Differentiation to appropriate fields thus yields the scaling behavior [4] as 

Ato(iM c ,T c ,L)=A®+gL* + ---. (13) 

Since the irrelevant exponent yi is in general not equal to — \2y t — d\, comparison with Eq. (12) 
shows that the last two terms in Eq. (12) must cancel one another, i.e., 

yi = -\2y t -d\. (14) 
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This accounts for the new correction exponent in Eq. (9). Further, if the asymptotic value 
is universal, one expects that this universal value applies to the canonical ensemble as 
well if 2y t -d< 0. 

We test these predictions numerically, by means of Monte Carlo simulations of five models 
with periodic boundary conditions, belonging to the d — 2, q — 1, 2 and 3 Potts universality 
classes. These models are: 

(1) the hard-square lattice gas model 

n = -Kj2^j-^J2 ak ( (T = 1 '°)' ( 15 ) 

(nn) k 

where the sum is on all nearest-neighbor pairs. The particles a = 1 have a 'hard' core so 
that nearest-neighbor exclusion applies, i.e., K — > — oo. The critical chemical potential is [5] 
p c = 1.334 015 100 277 74(1) with particle density p c = 0.367 742 999 041 0(3). 

(2) the Blume-Capel model [6] on the square lattice 

H=-KY,s i s j -»Y, s l (s = 0,±l). (16) 

(nn) k 

Zero spins s = are called 'vacancies'. For K — 1 the critical point [7] lies at p c (K = 1) = 
-1.702 7178(2), with a vacancy density p vc = 0.349 583 0(2). 

(3) the dilute 3-state Potts model on the square lattice 

H = -KJ2 &w (1 - k„o) + f E <W (17) 

(nn) i 

Each site carries a Potts variable a — 1, • • • , q, or a vacancy a = 0. The chemical potential 
p: controls the number of vacancies (N v ). A critical point was located [7] at K — 1.16940, 
fi c = -1.376 413(4), with a vacancy density p vc = 0.105 273 0(2). 

(4) Baxter's hard-hexagon model [8], also described by Eq. (15) with K — > — oo, but defined 
on the triangular lattice. Its critical point is [8] p c = ln[(ll + 5\/5)/2], with particle density 

Pc = (5 - VE)/io. 

(5) the bond percolation model on the square lattice, whose percolation threshold occurs at 
bond probability p c = 1/2. 
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For these models we sampled dimensionless quantities Q defined on the probability dis- 
tribution of the order parameter, which is the magnetization m for Isinglike models such as 
the Blume-Capel model. For that case we employ the Binder ratio [2] 



For the hard-square model, the same definition of Q applies with m replaced by the difference 
pi — p 2 of the particle densities on the two sublattices. For the dilute 3-state Potts model and 
the hard- hexagon model we used Eq. (18) with m 2 replaced by P1+P2+P3 — P1P2 — P2P3 — Pspi 
where pi denotes the density of the Potts variables in state i, and the particle density on 
sublattice % respectively. The mapping of the Ising model on the random-cluster model [9] 
makes it possible to express Q in moments of the cluster size distribution as 



with the moments of the cluster size distribution defined as l k = L~ 2k J2i=i c i where q is 
the size of the ith cluster, and N d is the total number of clusters. Eq. (19) defines Q for the 
bond percolation model. 

The conservation of particles in the canonical ensemble requires a definition of a 'particle'. 
This meaning is obvious for both lattice gases. In the percolation case, occupied bonds serve 
as 'particles': our simulations conserve the number of occupied bonds as N b = L 2 . This is 
simply realized by Monte Carlo moves that randomly interchange an occupied bond and an 
empty one. For the Blume-Capel and the Potts model we interpret nonzero spins Sj 7^ 
and variables <Tj 7^ as particles. All simulations, except for the bond percolation model, 
used the geometric Monte Carlo algorithm [10] that is well suited because it conserves the 
particle number and suppresses critical slowing down. For the Blume-Capel and the dilute 
Potts model it was combined with Wolff [11] cluster steps that do not affect the vacancies. 

For the bond percolation model, L 2 /A pairs of edges were interchanged between two 
subsequent samples. The sampling procedure decomposed the whole lattice into percolating 
clusters. The system sizes took 20 values in the range 4 < L < 4000. About 4 x 10 7 samples 



Q = ( m 2 ) 2 /(m 4 ). 



(18) 



Q = {h) 2 /wi)-m] 



(19) 
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were taken for each L for L < 800, and 4 x 10 6 samples for each L for L > 800. For the 
hard-square lattice gas, we simulated systems up to L = 4000, and for the hard-hexagon 
gas up to L = 1680. We took about 2 x 10 7 or more samples for each system size. For the 
Blume-Capel model and the q = 3 dilute Potts model, we used systems up to L = 1600 and 
L = 800, respectively; about 2 x 10 7 or more samples were taken for each system size. Since 
in a finite system the critical number of the particles p c L 2 is not an integer, we obtained the 
data at criticality by interpolation. For each of these models we determined Q at criticality, 
using only modest computer resources. 

Fig. 1 compares the finite-size corrections in Q of the canonical and the grand canonical 
hard-square lattice gas. The corrections are much more prominent in the canonical case, 
and decay more slowly with L. In the grand ensemble, Q scales as Q = Q c + biL Vi + b a L Va 
with an irrelevant Ising exponent yi = —2 and an exponent y a = 2 — 2y h = —7/4 describing 
the analytic background in (m 2 ). The canonical Q(L) data are shown versus 1/logL in 
Fig. 2; the approximate linearity indicates that the leading finite-size dependence is of a 
logarithmic nature. This is in line with Eq. (14), which yields y± = for Isinglike models. 
We thus fitted the Q data by 

2 

Q = Qc + $> fc (hiL - lnL )- fc + biL* + b a L y % (20) 

k=l 

where b^ and L are constants, with fixed exponents yi = —2 and y a = 2 — 2y h = —7/4. 
The Q data in the range 8 < L < 4000 are indeed well described by Eq. (20). We obtain 
Q c = 0.8563(2), L = 3.0(1), h = 0.190(2), and b 2 = -0.102(6); the value of Q c is in 
agreement with the grand canonical value Q c = 0.856216(1) [12,13]. We also fitted without 
logarithmic corrections, by 

Q(L) = Q C + ftiL* + 6 L»- + b 2 L^ + b 3 L^-\ (21) 

which was unsatisfactory because only a narrow range of sizes L could be accommodated, 
while the estimate y\ = —0.15(3) of the exponent is close to 0. Also for the Blume-Capel 
model we observed slowly decaying finite-size dependence in Q. The Q data for 12 < L < 
1600 were satisfactorily fitted by Eq. (20). We obtain Q c = 0.8567(7) and L = 3.3(1). 
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The Q data for the hard- hexagon model were fitted by Eq. (21) with exponents y a = 
2-2y h = -26/15 and y { = -4/5 [14]. The fits yield Q c = 0.9223(5) and Vl = -0.403(6) for 
the lattice gas. We have also determined the grand canonical value Q c = 0.85670(4), clearly 
different from the canonical value. For a similar fit to the data for the q = 3 dilute Potts 
model we find Q c = 0.9190(5) and Vl = -0.408(10). 

The Q data for the bond-percolation model were fitted by Eq. (21) with exponents y a = 
2-2y h = -43/24 and y t = -2 [14]. We obtain Q c = 0.87052(8) and y x = -0.499(2) « -1/2. 
For the grand canonical version of this model, we find [15] that the Q(L) data are also well 
described by Eq. (21), but without the term with exponent y± = —1/2. The universal ratio 
Q c = 0.87053(2) [15] agrees well with the canonical value. 

Our numerical analyses of Q, together with the predictions of Eq. (14), are summarized 
in Table I. 

Finite-size effects with exponent y± = —\2y t — d\ are not restricted to dimensionless 
quantities. Our derivation applies generally to the leading singular finite-size-dependent 
terms, which are modified by a similar factor as in the case of Q. We tested this numerically 
for the susceptibility (or its analog) of each of the the five models. Indeed we find effects in 
the critical canonical susceptibilities x^ according to x^{L) = x^(L)(l + rL yi + •■■). A 
further verification concerns the cluster- number density p c \ = N c i/N of the bond-percolation 

model, which at criticality behaves as [16] p c i(L) = p ca + L~ 2 (b + biL~ 2 H ) in the grand 

canonical ensemble, with p ca = (3v^3 — 5)/2 [16]. The p c i{L) data in the canonical ensemble, 
shown in Fig. 3, agree with y\ = —1/2. The fit of the canonical data by 

Pd(L) = p ca + L- 2 (b + hL^ + b 2 L~ 2 ), (22) 

confirms this value as y\ = —0.51(1) ~ —1/2. 

As another generalization, we expect finite-size effects with exponent — \2y h —d\ in systems 
with a conserved magnetization. 

The authors are indebted to Michael E. Fisher, Jouke R. Heringa, and Xiaofeng Qian for 
valuable discussions. This research was supported by the Dutch FOM foundation ("Stichting 
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FIGURES 
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FIG. 1. Finite-size dependence of Q of the canonical (A) and the grand canonical (□) 
hard-square lattice gas. 
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FIG. 2. Ratio Q of the canonical hard-square gas vs. 1/lnL. 
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FIG. 3. Scaled finite-size dependence of the cluster-number density p c \ of the canonical 
bond-percolation model vs. L -1 / 2 . 
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TABLES 

TABLE I. Summary of the main numerical results. The superscripts (n) and (t) indicate the 
numerical estimate and the theoretical prediction, respectively. 



model 


(n) 


(t) 




Qi 9) 


hard square 








0.8563 (2) 


0.856216 (1) 


Blume-Capel 


« 





0.8567 (8) 


0.856216 (1) 


q = 3 Potts 


-0.408 (10) 


-2/5 


0.9190 (5) 




hard hexagon 


-0.403 (6) 


-2/5 


0.9223 (5) 


0.85670 (4) 


bond percolation 


-0.499 (2) 


-1/2 


0.8705 (1) 


0.87053 (2) 
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